Star-shaped Local Density of States around Vortices in a Type II Superconductor 



Nobuhiko Hayashi,0 Masanori Ichioka,i and Kazushige MachidaS 
Department of Physics, Okayama University, 
Okayama 700, Japan 

The electronic structure of vortices in a type II superconductor is analyzed within the quasi-classical 
Eilenberger framework. The possible origin of a sixfold "star" shape of the local density of states, 
observed by scanning tunneling microscope experiments on NbSe2, is examined in the light of the 
three effects; the anisotropic pairing, the vortex lattice, and the anisotropic density of states at the 
Fermi surface. Outstanding features of split parallel rays of this star are well explained in terms 
of an anisotropic s-wave pairing. This reveals a rich internal electronic structure associated with a 
vortex core. 



On 

3 



(N 



> . 
(N ■ 

O ■ 

o : 
\o . 
o ■ 
\o • 
o\ : 

^ : 
o 
o . 

u '. 

3 ■ 



Much attention has been focused on a vortex struc- 
ture in high T c cuprate superconductors, yet some of 
fundamental problems associated with vortices still re- 
main even in conventional superconductors. The elec- 
tronic structure around a vortex core is one of them. No- 
tably, Hess et al. |[]-[| have done a series of remarkable 
and beautiful scanning tunneling microscope (STM) ex- 
periments on a layered hexagonal compound 2H-NbSe2 
(T C =7.3K) to reveal the detailed spatially resolved elec- 
tronic structure around a vortex core. They provide di- 
rect and vivid images of individual vortices and the flux 
line lattice. It was shown theoretically that the STM 
images observed by Hess et al. and Renner et al. Q re- 
flect the local density of states (LDOS) for bound states 
around a vortex. On the basis of the Bogoliubov-de 
Gennes (BdG) theory |q-[lO[| or the quasi-classical Eilen- 
berger (QCE) theory ]11| , |12| , the case of the cylindrical 
symmetric isolated vortex was considered there. How- 
ever, the detailed structure of the LDOS observed is not 
fully understood. 

The outstanding features of the experimental findings 
[§-§1 are summarized as follows when the magnetic field 
H is applied perpendicular to the hexagonal plane: (1) 
The LDOS for quasiparticle excitations has a sixfold 
"star" shape centered at a core. (2) The orientation of 
this star depends on the bias energy. At zero bias the 
"ray" extends away from the nearest-neighbor vortex di- 
rection [Fig. 0(a)] where the conventional 60° Abrikosov 
vortex lattice is formed. Upon increasing the bias voltage 
the star rotates by 30° [Fig. [|(c)]. (3) In the intermedi- 
ate bias voltage a ray splits into a pair of nearly parallel 
rays, keeping its direction fixed [Fig. [3(b)]. (4) In the 
spectral evolution as shown in Fig. 9 in Ref. |(| (also 
see Fig. 6 in Ref. ||) there exist additional inner ridges 
other than the outer ridges evolving from the zero bias 
peak at a core into the bulk BCS like gap edges far from 
a core. 

To understand these experimental results, Gygi and 
Schliiter jl(||l3| discussed the sixfold symmetric structure 
on the basis of the BdG theory by introducing a sixfold 
symmetric perturbation term and using their numerical 
solutions of the cylindrical symmetric case. While they 



explained the above mentioned features (1) and (2), the 
origin of their perturbation term is not clear, and it is 
uncertain whether the sixfold symmetric term can be 
treated by a perturbation theory. We notice that the 
features (3) and (4) are difficult to be explained in their 
framework. In the following, we show that an anisotropic 
pairing (or gap) effect which has been ignored so far is 
crucial for the bound states around a vortex. 

The purpose of this paper is to clarify the origin of 
the characteristic star-shaped LDOS by the QCE the- 
ory The QCE approach is more suitable than the 
BdG approach for the case cylindrical symmetry is bro- 
ken. The anisotropic pairing can be treated easily by the 
QCE theory. As for the origin of the sixfold symmet- 
ric structure, the following possibilities are enumerated; 
(A) the effect of an anisotropic pairing, (B) the effect of 
nearest-neighbor vortices, that is, the effect of the vortex 
lattice, and (C) the effect of the anisotropic density of 
states at the Fermi surface (FS). We calculate the LDOS 
for each case (A)-(C), and investigate which case can ex- 
plain the LDOS of the STM experiments. 

We introduce an anisotropic s-wave pairing A(r, k) = 
Ao(r)F(0) — Ao(r)(l + cacos66*) with hexagonal sym- 
metry and assume two-dimensional circular FS, neglect- 
ing a small warping along the hexagonal z axis (|| H), 
which is appropriate to NbSe2- Here r = (x,y) is the 
center of mass coordinate and k the relative coordinate 
of a Cooper pair. Now k is denoted by an angle 9 mea- 
sured from the a axis (or x axis) in the hexagonal plane. 
Thus the parameter ca denotes the degree of anisotropy 
in a pairing [ p^|Jl6| . The case ca = corresponds to a 
conventional isotropic pairing. 

Our calculation is performed after the method of Ref. 
for the isolated vortex and that of Klein ffj,(l9| for 
the vortex lattice. We consider a transport-like Eilcn- 
for the quasi-classical Green function 
9 if 



berger equation |l4j 
g{r, 9,iuj n ) = -in 
namely, 
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in a 2x2 matrix form, 
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supplemented by the normalization condi- 
tion: g(r,8,iu> n ) 2 = —ir 2 l. The Fermi velocity is vp. 
The bracket [ , ] is a commutator. The self-consistent 
equation is given by 

f r\9 t 

A(r,0) = 2^r^ ^- P (e')v(e,e')f(v,e',iu n ), (2) 

where the pairing interaction V(9, 9') = vF(9)F(9') is as- 
sumed to be separable. The ^-dependent density of states 
at the FS is p(9). We calculate A (r) self-consistently 
from Eqs. (1) and (2) for T/T c = 0.1. The LDOS is 
given by 

N(r,E) = J — p (9)N(r,9,E) 

= J —p(9)Reg(r,9,iu Jn ^E + ir 1 ), (3) 

where r\ (>0) is an infinitesimal constant. To obtain 
g(r, 9, iuj n — > E + irj), we solve Eq. (1) for rj — iE instead 
of the Matsubara frequency ui n using the self-consistently 
obtained pair potential. From now on, the energy E and 
the length r are scaled by the uniform gap Ao at zero 
temperature and £ = vf / A respectively. 

Let us start out by considering (A) the effect of an 
anisotropic pairing. The anisotropic pairing is suggested 
from the STM experiment at zero field. The observed 
I-V tunneling spectrum in the bulk indicates a substan- 
tial gap anisotropy (the gap amplitude with the averaged 
value l.lmV distributes from 0.7mV to 1.4mV, see Fig.l 
in Ref. ||), which is consistent with the density of states 
in the anisotropic pairing case (see Fig. 5 in Ref. |l6|]). 
To see the influence of (A) clearly, we neglect the other 
effects (B) and (C) here. We calculate the LDOS in the 
isolated vortex case assuming an isotropic FS. 

In Figs. 1(d), 1(e), and 1(f), we show the calculated 
LDOS for several bias energies E in the case ca = 0.5. 
It is seen from Fig. 0(d) that the sixfold star centered at 
a core is oriented away from the x axis by 30° for E=0. 
The rays along the 30° direction can be seen to extend up 
to r ~ 3 from the center. It corresponds to the following 
fact. As for the bound state around a vortex core at 
E = 0, the wave function of the quasiparticle running 
through the vortex center has the maximum amplitude. 
Among them, because F(9) is most depressed at 9 = 30°, 
the quasiparticle with the momentum k along the 30° 
direction feels weakest superconducting pair potential, 
and its wave function keeps large amplitude until far from 
the vortex. 

It is seen from Fig. 0(e) that at the intermediate energy 
each ray splits into two parallel rays, keeping its direc- 
tion. This characteristic feature is precisely observed by 
Hess [Fig. 0(b)] (also see Ref s. [flp|] ) and is absent in the 
result by Gygi and Schliiter Jl0| , |13| . For this intermedi- 
ate excitation energy, the wave function of the quasipar- 
ticle passing at a small finite distance from the vortex 



center has the maximum amplitude. Among them, the 
quasiparticle with the momentum k along the 30° and 
its equivalent directions keeps large amplitude until far 
from the vortex, and forms the nearly parallel structure 
in the r space. 

With increasing the bias energy E further, the sixfold 
star becomes a more extended one, and its orientation 
rotates by 30° as seen from Fig. 0(f). Note that the heads 
of each ray spread out. It coincides with the observation 
in Fig. 0(c). 

Another way to examine the quasiparticle excitations 
in the vortex state is to see how the spectrum evolves 
along the radial lines. The calculated spectrum evolution 
shown in Fig. ^ is compared with the data in Fig. 9 in 
Ref. [|| (also see Fig. 6 in Ref. ||). Here the radial line is 
chosen to make an angle 15° from the x axis, where there 
are five ridges for E > 0. We label these ridges as as. As 
r increases, two higher energy ridges a and (3 approach 
the maximum of the gap energy 1 + ca, and the ridge 
7 is found to become very weak intensity. The ridge 5, 
which approaches the minimum of the gap energy 1 — ca 
with increasing r, corresponds to that reported in the 
isotropic pairing model by Gygi and Schliiter (see Fig. 15 
in Ref. JkJ). The ridge e is the result reflecting the sixfold 
symmetric ray structure. The ridges 8 and e correspond 
to the outer and inner ridges in the experimental data 
HH by Hess et at, respectively. 

To clarify the behavior of these ridges, we show in Fig. 
U the trace of the ridges in the LDOS as a function of r 
along the radial lines for 0° (a), 15° (b), and 30° (c) from 
the x axis. For comparison, the experimental data || are 
also presented. We direct our attention to the inner ridge 
(line e), which strongly depends on the radial angle. It 
moves toward the lower energy side as the angle increases. 
At 30°, it reduces to the ridge at E = (horizontal r- 
axis). This behavior of the inner ridge agrees well with 
that of the experimental data. 

Thus, the effect of the anisotropic pairing can explain 
not only (1) the sixfold star shape and (2) the 30° rota- 
tion, but also (3) the split parallel ray structure at the 
intermediate energy and (4) the characteristic behavior 
of the inner ridges in the spectral evolution. 

Notice that this characteristic spreading of the ray 
structure depends on the degree of anisotropy ca; In the 
isotropic limit (ca = 0) the star shape turns into a cir- 
cular object. As ca increases, the star structure spreads 
out farther from a core. This implies that in NbSe2 a 
substantial gap anisotropy must exist even within the 
hexagonal plane. It is consistent with the tunneling ex- 
periment || in the bulk mentioned earlier. 

Let us argue a possible origin of this rather strong gap 
anisotropy. It is known that the system enters into a 
charge density wave (CDW) state at 32K which coexists 
with superconductivity at lower temperatures. The FS's 
consist of two kinds of cylinders, each centered at the F 
point and K point in reciprocal space pfl] . The CDW re- 
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organizes the FS's, opening up the CDW gaps by nesting. 
The remaining FS's left over by the CDW gap formation 
are available for the superconductivity. Therefore, we can 
expect a substantial gap anisotropy within the hexagonal 
plane. 

Next, we turn to (B) the effect of the vortex lattice. 
The LDOS is calculated in the case of a triangular vortex 
lattice by using the self-consistently obtained pair poten- 
tial. Here, neglecting the other effects (A) and (C), we 
consider the isotropic s-wave pairing (ca = 0) and as- 
sume an isotropic FS. As for the LDOS in the vortex lat- 
tice case, Klein calculated only the momentum-resolved 
LDOS N(r,8,E) in Eq. (3) for specific k direction par- 
allel to mri + riY2 (m and n are integer, ri and V2 are 
unit vectors of the vortex lattice) because the so-called 
symmetric method was used fll| . Using the explosion 
method |l7|-[l9[|, we succeed to calculate N(r,0,E) for 
arbitrary k direction. Then we can obtain the LDOS 
integrated over k directions. 

Important results obtained for the vortex lattice are; 

(i) We do find a sixfold star shape of the LDOS and the 
30° rotation upon elevating the bias energy. The sixfold 
star originates from the triangular vortex lattice effect. 

(ii) The orientation of the star coincides with the STM 
data, namely the ray extends toward the next nearest- 
neighbor vortex direction at lower bias energy. Therefore 
we succeed in determining the absolute direction relative 
to the flux lattice. This is one of the most eminent fea- 
tures in the STM data and absent in Gygi and Schliiter 
Jl0|,[l3]]. (iii) However, the split parallel ray structure at 
the intermediate energy is not reproduced, (iv) The char- 
acteristic sixfold symmetric LDOS appears only at a high 
magnetic field such as 1 Tesla for the material parame- 
ters appropriate to NbSe2 (The BCS coherence length 
77 A and the BCS penetration depth 690A where 
the core regions substantially overlap each other. At a 
lower magnetic field such as 0.1 Tesla, the LDOS reduces 
to the almost circular structure. 

As for (C) the effect of the anisotropic density of 
states at the FS, we find that reflecting the assumed six- 
fold anisotropy the LDOS exhibits the characteristic star 
shape and the 30° rotation even for isotropic s-pairing. 
However, the split parallel ray structure at the interme- 
diate energy is not reproduced. 

In conclusion, by solving the quasi-classical Eilenberger 
equation self-consistently, we have examined the possible 
origins of the star-shaped LDOS observed in the STM ex- 
periments on NbSe2. Among the three effects (A), (B), 
and (C), each one can reproduce the experimental fea- 
tures (1) the sixfold star shape and (2) the 30° rotation. 
In NbSe2, all the three effects may play important roles. 
However, to reproduce the features (3) the split paral- 
lel ray structure at the intermediate energy and (4) the 
characteristic behavior of the inner ridges in the spectral 
evolution, we have to consider the LDOS by including 
the effect of the anisotropic pairing. 



As for the relation between the orientation of the star 
shape and that of the vortex lattice, the calculation 
for the vortex lattice has to be done by including the 
anisotropy effect. It is noted that the orientations of the 
star shape and of the vortex lattice relative to the un- 
derlying crystal lattice are internally correlated through 
the gap function. The gap function itself is determined 
by the FS topology. The complete answer to this can 
be obtained by considering all the three effects simulta- 
neously, which belongs to a future work. The present 
study underlies the basic understanding for fundamental 
electronic properties of vortices both for conventional and 
unconventional superconductors, revealing a rich internal 
electronic structure associated with vortices. 

In connection to the anisotropic pairing effect, we point 
out that high T c superconductors should exhibit similar 
STM characteristics if a <i-wave pairing such as A(r, k) = 
Ao (r) cos 29 is realized. We predict a "fourfold" star and 
its 45° rotation there 0. 
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present project. 
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FIG. 1. Tunneling conductance images observed by Hess 
et al. at 0.1 Tesla for the bias voltage 0.0m V (a), 0.24mV (b), 
0.48mV (c), where 1759A x 1759A is shown (also see Refs. 
[4,5]). The nearest-neighbor direction of the vortex lattice 
is the horizontal direction. The LDOS images calculated for 
E = (d), 0.2 (e), and 0.32 (f), where 6£ x 6£ is shown. 

FIG. 2. Evolution of the spectra along a radial line at the 
angle 15° from the x axis. Each ridge is labeled as as from 
the high energy side. 

FIG. 3. The ridge energies of the LDOS as a function of r 
along the radial lines for 0° (a), 15° (b), and 30° (c) from the 
x axis (solid lines). The case (b) corresponds to Fig. 2, where 
the labels a-e are the same as those in Fig. 2. Experimental 
data[3] are also presented by points • (outer) and o (inner), 
where r and E are scaled by 350A and 1.67mV, respectively. 
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